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Abstract 



The resonant part of the B\ g electronic Raman scattering response is cal- 
culated within the t — J model on a planar lattice as a function of temperature 
and hole doping, using a finite-temperature diagonalization method for small 
systems. Results, directly applicable to experiments on cuprates, reveal on 
doping a very pronounced increase of the width of the two-magnon Raman 
peak, accompanied by a decrease of the total intensity. At the same time the 

peak position does not shift substantially in the underdoped regime. 
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Proper understanding of electronic properties of cuprates, representing a remarkable 
example of strongly correlated systems, remains a challenge for theoreticians and experi- 
mentalists. One of very useful probes for the investigation of electronic excitations in these 
materials has been the Raman scattering. Using the latter method, it has been clearly es- 
tablished that in the reference insulating substance, e.g. in La 2 Cu0 4 , the most pronounced 
electronic Raman processes at low temperatures correspond to the short-wavelength mag- 
netic (two-magnon) excitations Q, which can be well described within the Heisenberg model 
for the planar antiferromagnet (AFM) 0. 

A more general framework for the theoretical explanation of the Raman scattering in 
strongly correlated systems, and specifically in cuprates, has been first given within the 
Hubbard model, where the effective Raman operator for resonant and off-resonant condi- 
tions has been derived || . Nevertheless, due to the difficulties in performing the theoretical 
analysis and due to the possibility of other relevant processes, an agreement on the appro- 
priate interpretation is still lacking for a number of Raman scattering results in cuprates. 

One aspect concerns the pronounced temperature dependence of the linewidth in the 
undoped AFM The latter has been attributed to the phonon-induced broadening ||, but 
an important role could be as well played by higher-order resonant processes 0, becoming 
relevant in the resonant-Raman conditions. 

Another problem is the doping dependence on the Raman scattering |7||§. Recent ex- 
periments, performed on YBaCuO materials in the resonant regime ||, show a dramatic 
increase of the broadening of the two-magnon peak with doping, so that spectra appear 
essentially flat in the normal phase T > T c when approaching the 'optimum' doping cases 
with highest T c . At the same time, the peak position does not move appreciably. 

The aim of the present paper is to analyse the influence of doping and finite temperatures 
on the Raman spectra within the framework of the one-band t — J model, assuming that the 
dominant resonant contribution remains of the Loudon-Fleury type due to the spin exchange 
0. Assuming the Raman operator we calculate the spectra at various hole concentrations 
and temperatures T by employing the recently introduced finite-T diagonalization method 
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for small correlated systems ||10|| . T > results are of interest by themselves. However, in 
the case of Raman spectra, where quite broad features are expected particularly in doped 
systems, using finite T > mainly represents a technical advantage to obtain macroscopic- 
like spectra even for small systems, as already established in previous applications of the 
method to the evaluation of optical conductivity and dynamical spin correlations [JCT]]. The 
improvement over the usual T = Lanczos method can be judged by comparing our results 
with previous Raman spectra obtained within the t — J model containing few mobile holes 



E| . While some aspects appear similar (e.g. lowest frequency moments), actual spectra 
calculated at T = from the ground state are dominated by few peaks, which are clearly 
size-dependent. 

One of the most studied prototype models of correlated systems, and specifically of the 



low-energy properties of cuprates, is the t — J model [13 



H = -t^2 ( c )s c is + H.c.) + J XX<% ■ Sj - -UiUj), (1) 

(ij) 

where cj s (cj S ) are projected fermionic operators, prohibiting double occupancy of sites. The 
ground state of the model has been studied both by analytical and numerical methods 



rj| . Recent investigations of the finite-T dynamical spin and charge response [[D]] established 



also qualitative (as well as quantitative) contact with anomalous normal-state properties of 
cuprates, confirming the t — J model as the proper framework for studies of the low-energy 
physics for such doped AFM. 

The operator, relevant for the resonant Raman scattering in cuprates, cannot be deter- 
mined uniquely within the t — J model (in contrast to operators for some other charge and 
spin response functions), since necessarily Raman processes involve higher resonant levels. 
Here we adopt a view that more complete models for cuprates, e.g. the three-band model, 



can be for processes of interest mapped onto an effective t — U Hubbard model [|15| , |I"3"f , where 
the upper Hubbard-band states now correspond to charge transfer excitations in cuprates. 
Within the Hubbard model near half-filling the Raman-scattering operator has been derived 
in the limit t/U < 1 [[|, yielding the well known form for the Heisenberg AFM || 



Ri = A^2(e inc ■ fij)(e sc ■ rij)(Si ■ Sj - ju^), (2) 

(ij) 

where e inc , e sc are the incident and the scattered electric-field vector directions, respectively, 
fij is the vector connecting sites i and j, and the amplitude factor A = At 2 /{U — u>i nc ) oc J 
incorporates the resonance at the incident-light frequency uji nc ~ U. 

For the undoped (insulating) case the operator Ri is the only resonant term of the 
order t 2 /U, hence the dominant one, at least outside the resonance \U — uj inc \ > t. In the 
doped AFM additional processes of the same order of magnitude are possible. The latter 
correspond to the hopping term involving three neighboring sites, leading partly also to the 
next-neighbor hopping correction, i.e. to the t' term, within the t — t' — J model [|I3| ,|15|1, 

(i ,j,k),ss' 



R2 — J]] ( e inc ' r ij){tsc ' r jk) c ks' C js /C Js C is- (3) 



The contribution of such term clearly scales with the hole concentration (it vanishes at half 
filling), so it is a quantitative question whether it can become important in the relevant 
parameter regime of cuprates. 

The Raman spectral function is given by 



1 f°° 

Uu) = Re / dt e luJt (R(t)R(0)) 

rrN Jo 



1 J2e-^\(m\R\n)\ 2 5(u-E m + E n ) , (4) 



NZ 

- 1 * n,m 

where u = u sc — uj inc , (3 = 1/T (we use furtheron ks = h = 1), and Z is the partition 
function. Sums run over all eigenstates n, m with corresponding energies E n , E m and are 
size-normalized with 1/N. 

To calculate I{uS) within the t — J model on a square lattice we study finite clusters 
with N = 16 — 20 sites. Via Eq.(P|) this would be impossible with the method of full 



diagonalization ||Iq| , used previously to investigate the Raman response for the Heisenberg 
model with maximum lattice size N = 16. Instead, we employ the recently developed method 
for finite-T dynamical (and statical) correlation functions [|10|| , based on the Lanczos iteration 



procedure combined with random sampling. The method is quite convenient for the analysis 



of the Raman response, in particular for doped systems, since quite broad spectral features 
and a weak T-dependence require a modest number of Lanczos steps M < 100 ||iO| , |rT|| . In 
results presented below we thus reach systems with N st ~ 400.000 states in the largest q, S z - 
basis sector, while the sampling involves typically Nq ~ 200 — 500 initial configurations. As 
discussed quite extensively in connection with previous applications, the spectra reveal a 
macroscopic-like behavior only at finite T > T*, where T* is related to the low-energy level 
spacing and thus dependent on the system size and hole doping. In all calculations we will 



fix J/t = 0.3, as relevant for cuprates [13|, where J ~ 950 cm -1 . For such parameters we 
typically reach T* ~ 0.1 t = J/3 in the doped system with 0.1 < Ch < 0.3. 

We restrict our analysis to the dominant B Xg scattering geometry with e inc = (e x +e y ) / \[2 
and e sc = (e x — e y )/\/2. Let us first discuss results for the undoped Heisenberg model. Here 
the low-energy levels are quite sparse and hence we reach only T* ~ 0.5 J in the largest 
N = 20 system. In Fig. 1 we present the spectra for several T > T*. Consistent with 
T ~ results |2||J and with the low-T experiments we observe a two-magnon peak at 
uj ~ 3.3 J, being quite narrow at low T < J, and having the width limited by quantum 
fluctuations 0. From Fig. 1 it follows that the peak width starts to increase substantially 
only at higher T ~ J, where a gradual transition to a broad featureless spectrum occurs. 
It should be also noted that in spite of large broadening the peak does not move with T 
up to T ~ 2 J. Nevertheless, when considering undoped cuprates, such as La 2 Cu0 4 with 
J ~ 1400 K, clearly other mechanisms have to be invoked to account for the observed 
pronounced T-dependent width at lower T ^ J [||||. 

We proceed to the doped systems, where we first consider only the exchange part R = Ri, 
Eq.(0). Here T seems to play a less essential role, provided that T* < T < J. Typically 
we observe for doped systems only a steady decrease of the intensity with T, e.g. a ~ 20% 
reduction in the interval T/J = 0.3 — 1.0 for Ch ~ 0.2, and ~ 30% for Ch ~ 0.05. On the other 
hand, the dependence on doping is essential, as evident from Fig. 2, where we present spectra 
for various hole concentrations Ch = N^/N < 0.25 at lowest T = 0.15 t > T*. Already the 
smallest possible (for available sizes N) finite doping Ch = 0.05 increases dramatically the 



width of the two-magnon peak, and spectral features become overdamped on approaching 
the 'optimum' doping c/j ~ 0.2. It is however remarkable, that the peak position does not 
shift appreciably in the underdoped regime Ch < 0.1. Only for Ch > 0.15 the spectra change 
to a broad central-peak form with a maximum at uj = 0. 

Another experimentally relevant quantity is the total spectral intensity Jo = / °° I(u)du 
and its variation with doping. At the same time, we evaluate the frequency moments (u), 
(uj 2 ), and the spectral width o = <J {uj 2 ) — (uj) 2 , where the averages are calculated with 
respect to the I(ui > 0) part of the spectra. Quantities on Fig. 3 are calculated within 
various-size systems for T = 0.5 J, where values should be quite close to the T = results. 
Results for finite > are presented for N = 16 — 20. For undoped system N h = we 
include values for N = 20, as well as for N = 26 at T = 0. In the latter case, we get e.g. the 
width a ~ 0.8 J, very close to analytical predictions 0. The scattering of values in Fig. 3 
for similar Ch could be partly attributed to the error due to the restricted random sampling 
pTjj] for larger systems, while smaller systems at chosen T ~ T* could still suffer from finite 
size effects. 

As evident from Fig. 3, the total intensity J decreases steadily but substantially with 
doping, i.e. by a factor ~ 3 from Ch = to Ch = 0.25. The average frequency (uj) starts to 
increase slowly only at Ch > 0.15. The enhancement can be attributed to the emergence of 
high-frequency tails with u > t in doped systems, whereas the peak position moves in the 
opposite way. On the other hand, the spectral width a shows a more dramatic variation, in 
particular at low doping < 0.1. 

Let us finally discuss the influence of the three-site hopping term R 2 , Eq.(^|). Our results 
indicate that it is less important at low T and in the concentration range of interest Ch < 0.25, 
either calculated separately as R = R 2 or combined R = R\+R 2 - When evaluated as R = R 2 
we observe that AI(u) is essentially featureless up to uj ~ 4 t, and becoming larger only at 
T > J, while at T ~ T* it amounts only to e.g. 0.2 Jo at Ch = 2/16. Taken into account 
as R = R\ + R 2 , the largest corrections are at Ch = 0.25 as expected, i.e. SIq/Iq ~ 0.3. 
Nevertheless, the most pronounced changes are not in the regime uj < 2 t (spectra are even 
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flatter there), but rather in the large-cu tails, leading to substantially increased (u>) (e.g. 
~ 70% for Ch ~ 0.25) and a at the 'optimum' doping. It seems that the second-neighbor- 
hopping processes are quite restricted at low T, due to remaining short-range magnetic 
correlations, disappearing only at higher T > J. It should be however noted that we did not 
introduce an analogous three-site term (or a simpler t' hopping term) in the t — J model, 
Eq.(|l]) |T^JT5|] , which possibly could lead to some quantitative modifications. 

How can we interpret the above results for the B\ g Raman scattering? Up to the over- 
doped situation Ch ~ 0.3 the scattering seems to be dominated by the spin-exchange part, 
Eq.(||), since the latter appears to determine the low-energy fluctuations even at the 'opti- 
mum' doping [11]. Still, there are evident changes with doping. In an undoped AFM a well 



developed longer-range order (long-range at T = 0) induces a sharp peak in In the 

underdoped regime Ch < 0.1 the reduced AFM correlation length does not shift the peak, 
but induces a large broadening and a reduction of the intensity. In this respect the doped 
system behaves quite similarly to an undoped AFM but at elevated T ~ J, as one can 
conclude from the similarity of spectra on Figs. 1,2. For Ch > 0.1 the spin system seems to 
approach a totally incoherent one. Any coherence in the contribution of different bonds in 
Eq.(^) is eliminated, leading to a reduction of the intensity and a broad featureless spectrum 
with a central maximum at u = 0. 

Finally let us comment on the relation of our results to experiments. A system- 
atic resonant-Raman scattering study has been recently performed on a sequence of 
YBa 2 Cu306+ x materials and YB^C^Og ||. Although it is not straightforward to relate 
in this case the actual doping of Cu0 2 layers with our model doping Ch (and consider also 
the possible role of CuO chains), it seems that experiments for the underdoped materials 
correspond well to our model results, both regarding the shape of the Raman spectra and 
their intensity variation with doping. On the other hand, entering the 'optimum'-doping 
and the overdoped regime both experiments and the theory show consistently a nearly flat 
Raman response. 

One of the authors (P.P.) wants to acknowledge stimulating discussions with M. Klein 



and G. Blumberg. 
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FIGURES 

FIG. 1. Raman intensity I /A 2 vs. uo/J for the undoped Heisenberg model at various T, as 
calculated for N = 20. An additional smoothening width A = 0.2 J is used at low T. 

FIG. 2. I /A 2 vs. lo/J for the t — J model at various hole concentrations and T = 0.5 J. 
The smoothening is A = 0.2 J. 

FIG. 3. Total Raman intensity Iq (triangles), the average frequency {w)/J (circles), and the 
spectral width a/ J (squares) vs. doping at T = 0.15 t, as evaluated for systems with N = 16—20. 
For Ch = also N = 26 results at T = are included. Lines are guides to the eye only. 
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